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ABSTRACT 


This study presents many aspects of contact 
problems and presents a solution for the determi- 
nation of the stresses and displacements in an arbi- 
trarily thick layer that is resting on an elastic 
(Winkler type) foundation. The contact problems dis- 
cussed are the nonlinear bending of beams, the stress 
in a thin plate resting on an elastic foundation, the 
stress in a half plane because of surface punches and 
pressures, the stress in a half space because of sur- 
face punches and pressures, and the stress in a layer 
resting on a rigid base because of surface punches 
and pressures. The problem of determining the 
stresses of a thick layer resting on an elastic founda- 
tion is then solved for a general surface pressure 
distribution. Numerical results are presented for 
the particular case of a surface pressure uniformly 
distributed over a circle. 


ii 



CONTENTS 


Section Page 

SUMMARY . . . . . 1 

INTRODUCTION 1 

SYMBOLS . . 2 

HISTORICAL BACKGROUND 6 

ELEMENTARY CONTACT PROBLEMS 7 

Nonlinear Bending of a Cantilever Beam 8 

Nonlinear Bending of a Simply Supported Beam 8 

Thin Plate Resting on an Elastic Foundation 8 

TWO-DIMENSIONAL PROBLEMS OF A HALF PLANE 9 

First Fundamental Problem for Half Plane . 11 

Punch Problems for the Half Plane . . 13 

FIRST FUNDAMENTAL AND PUNCH PROBLEMS FOR THE HALF SPACE ... 15 

First Fundamental Problem for a Half Space 17 

Punch Problem for the Half Space 20 

FIRST FUNDAMENTAL AND PUNCH PROBLEMS FOR AN ELASTIC LAYER 

RESTING ON A RIGID FOUNDATION 22 

First Fundamental Problem for an Elastic Layer Resting on a Rigid 

Foundation .................................. 23 

Punch Problem for an Elastic Layer Resting on a Rigid Foundation ..... 25 

SOLUTION FOR THE PROBLEM OF A THICK LAYER RESTING ON AN 

ELASTIC (WINKLER TYPE) FOUNDATION 26 

NUMERICAL RESULTS ON THE CALCULATION OF CONTACT STRESS 

BETWEEN A LAYER AND ITS ELASTIC FOUNDATION 36 

Derivation of the Contact Equation 38 

Discussion of Analytical Results ........................ 39 



Section 


Page 

CONCLUSION 40 

REFERENCES 42 


iv 



FIGURES 


Figure Page 

1 Nonlinear bending of a cantilever beam • • 44 

2 Nonlinear bending of a simply supported beam ............. 44 

3 Thin plate resting on an elastic foundation 45 

\ \ ... ; 

4 Pressure on bottom of thin plate resting on an elastic foundation 

and subjected to a concentrated load P .......... 45 

5 First fundamental problem 46 

6 Punch problem . . . 46 

7 Half plane loaded by a surface pressure 47 

8 Half plane loaded by a concentrated normal force 47 

9 Contact pressure between half plane and flat bottom punch . . . .... 48 

10 Contact pressure between half plane and a cylinder . . . . . . . . ... 49 

11 Contact pressure between a half planb and a wedge-shaped punch . . . 50 

12 Half space loaded by a surface pressure ........ . . . ..... 50 

13 Half space loaded by a concentrated normal force ........... 51 

14 Half space loaded by a pressure uniformly distributed over 

a circle 51 

15 Half space deformed by a punch . ...... 52 

16 Contact pressure between a half space and a flat-base cylindrical 

punch 53 

17 Contact pressure between a half space and a conical-base punch .... 54 

18 Elastic layer resting on rigid foundation and subjected to a 

surface pressure ............................ 55 

19 Elastic layer resting on rigid foundation and subjected to a 

pressure uniformly distributed over a circle 55 

20 Contact pressure between layer and rigid foundation (Sneddon’s 

solution for a layer on a rigid foundation) 56 


v 



Figure Page 

21 Elastic layer resting on rigid foundation deformed by a flat-ended 

cylindrical punch 56 

22 Thick layer on elastic foundation 

(a) Layer on elastic foundation subject to an axially symmetrical 


pressure 57 

(b) Cross section through axis of symmetry showing coordinate 

system ............................... 58 


23 Contact stress between a layer and its elastic foundation K ; the 
surface loading is a uniform pressure distributed over a circle 
of radius a 


(a) 

^ _ 0. 0001^ 

59 

2h(l - v) '• 


(b) 

■v _ 0. 001m 

. 60 

“ 2h(l - v) 


(c) 

K _ o. 01 m 

61 

2h(l - v) 


(d) 

rr _ 0. 1m 

62 

2h(l - v) 


(e) 

If _ U 

63 

2h(l - v) 


(f) 

Y - 

64 

2h(l - v) 


(g) 

TC — 

65 

2h(l - v) ; 


(h) 

y _ 10m 

. . 66 

2h(l - v) : 


(i) 

k= 1™E-. 

67 

2h(l - v) 


(j) 

v _ 1000m 

68 

2h(l - v) 


(k) 

K = 00 

69 



CONTACT PROBLEMS IN ELASTICITY AND THE SOLUTION FOR A 


THICK LAYER ON A WINKLER TYPE FOUNDATION* 

By William C. Schneider 
Manned Spacecraft Center 


SUMMARY 


This study presents many aspects of contact problems and presents the solution 
to the problem of an arbitrarily loaded thick layer resting on an elastic (Winkler type) 
foundation. When solving stress problems, the forces are usually transmitted by con- 
tact between two bodies. If the method of contact influences the stress, the problem is 
called a contact problem. Here, such problems are discussed in order of difficulty, 
beginning with the simplest. The first type of contact problem considered is the non- 
linear bending of beams. The solutions to the two-dimensional and three-dimensional 
problems are presented. Finally, the stress in a layer that is resting on a rigid foun- 
dation is determined. 

The problem of a thick layer resting on an elastic foundation is then solved to de- 
termine the effect of plate thickness on the stress distribution. The solution is given in 
terms of infinite integrals that are virtually impossible to integrate by hand. To plot 
the solution, the integrals are integrated numerically on a computer. Numerous plots 
are given for a particular surface loading which consists of a pressure that is uniformly 
distributed over a circle. 


INTRODUCTION 


This study has two main purposes. The first purpose is to present numerous 
examples of the many aspects of contact problems in elasticity. The second and more 
important purpose is to solve the problem of a thick layer resting on an elastic 
(Winkler type) foundation, the surface of which is subjected to any prescribed normal 
loading. Contact problems are of practical importance in the design of building founda- 
tions, ball bearings, clamps, or any machine part which involves forced contact be- 
tween two objects. The first section of the paper gives a brief history of the proposed 
solutions for the plane contact problems. The following four sections deal with aspects 
of punch problems; the solutions are presented, but the derivations are omitted. How- 
ever, the derivation of the solution to the problem of a thick plate resting on an elastic 
foundation is presented in detail in a later section because this solution is considered 
to be the main contribution of the paper. 


*This paper was originally written in partial fulfillment of the requirements for 
the master of science degree in mechanical engineering at the University of Houston. 



The advice and guidance of Dr. Z. O. Olesiak, Visiting Professor of Mechanical 
Engineering, University of Houston, Houston, Texas, is acknowledged.- 

SYMBOLS 


A arbitrary constants of integration associated with u^ 

A = J_ (n/2 - 1/2)1 
n " 'S* (l - /)(n/2 - 1)'. n 

a boundary of punch 

B arbitrary constant of integration associated with Wq 

b thickness of elastic layer resting on rigid foundation 

c distance from the origin to the force 

D plate rigidity 

d layer thickness 

E modulus of elasticity 

F total resultant force of pressure exerted on a surface 


’ (m, T)= J 


F 0 (m, t) = — 
2' ’ ' m 


w cosh wJQ(mw)Jj(rw) 


2w + sinh(2w) 


dw 


f 


cosh wJ 1 (mw)J 1 (Tw) 


dw 


2w + sinh 2w 

f^tp^ displacement under the punch for two-dimensional problems 
f^tgj constant part of deflection 

f^t^ varying part of deflection 


l^(e) Hankel transform of f(r) of the order v 
G = 2hK/l - v) 


2 



sinh wJQ(mw)Jj(Tw) 
2w + sinh 2w 


GjOn, t) 


G 2 (m, t) 


g(t) 

H,(» 


h 

1 

j 

v 

K 

K 0 

k 

kei 

ker 

L 

1 

M 

m 

P 

P 0 

P 0 

P 

p ( t o) 



^ j sinh wJj(mw)Jj(rw) 
m Jq 2w + sinh 2w(w) 

representation of solution of the Fredholm equation 
ith order Hankel transform 

half of plate thickness 

cross section area moment of inertia 

integrals of Lipschitz- Hankel type involving products of Bessel functions 
Bessel function of the first kind, order v 

modulus of foundation 

Bessel function of the second kind 

dummy variable 
Kelvin function 
Kelvin functions 
length of beam 
characteristic length 
bending moment on beam 
dummy variable 
concentrated normal load 
total vertical load 

zeroth order Hankel transform of the surface pressure 
uniformly distributed pressure 

normal load as a function of boundary horizontal coordinate 
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q 


normal load distribution on the surface 


R radius of curvature of rigid body 

Re real part of any function 

r radial (polar) coordinate 

T^t^ shear loading on the boundary 

t dummy variable 

tg x- coordinate of any point on the surface of an elastic body 

u displacement in x-direction (Cartesian coordinates) or in r-direction (polar 

coordinates) 

v displacement in y-direction 

w displacement in z-direction 

x, y, z axes of Cartesian coordinate system 

a angle between the vertical and the point of contact of a beam and a rigid 

body 

/? dummy variable 

y . dummy variable 

y shear strain on x-face in y-direction 

xy J 

5 beam displacement in the vertical direction 

6 deflection at point B (fig. 1) 

6^ Dirac delta function 

deflection of segment DB considered as simple cantilever beam (fig. 1) 

deflection at point B that results from the slope at point D (fig. 1) 

5_ deflection at point B representing the vertical distance of point D from the 

horizontal tangent at point A (fig. 1) 

e deflection of a beam 
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e x, y, z 

v 

e 


K 

X 

M 

v 

p 

CT 

r 

a 

rz 



a 

x 

a 

xy 

a 

y 

a 

z 



T 

X(t) 

* 


normal strain in x-, y-, or z-direction 

dummy variable equal to 2£h 
angular (cylindrical) coordinate 
constant equal to 3 - Au 
Lame constant 

Lame constant (shear modulus) 

Poisson ratio 

dummy variable 

radial (cylindrical) coordinate 

normal stress on radial face in radial direction 

shear stress on radial face in z-direction 

shear stress on radial face in 0 -direction 

normal stress on x-face in x-direction 

shear stress on x-face in y-direction 

normal stress on y-face in y-direction 

normal stress on z-face in z-direction 

normal stress on 0 -face in 0-direction 

normal stress on 0-face in ^-direction 

dummy variable 
complex function of z 
operational function of t 
dummy variable 
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operational function of £ 
any argument of H 


TO 

w 


Operators: 

( ) complex conjugate of a function 

(' ) derivative 

2 

v Laplaeian operator 


HISTORICAL BACKGROUND 


Sadovskii (ref. 1) solved several problems dealing with the pressure of a rigid 
body on an elastic semi- infinite plane. He considered the pressure of a plane-base 
punch on an elastic semiplane, and he also considered the case of infinitely many 
punches. The general method of solving plane contact problems by the Cauchy in- 
tegrals was given by Muskhelishvili in reference 2. Here, the case was established 
that the required function of a complex variable is regular everywhere, including the 
edges of the contact area. (Muskhelishvili' s Cauchy integral solutions are presented 
in a later section. ) Begiashvili (ref. 3) generalized the problem of one punch to solve 
a problem involving any number of punches. Lomidze (ref. 4) solved the problem con- 
cerning the pressure of a system of jointed ground foundations. Klubin (ref. 5) deter- 
mined the stresses arising inside a heavy elastic semiplane. A number of frictionless 
punch problems for a semiplane were considered by Shtaerman in reference 6. 

Numerous problems of elastic equilibrium, both two-dimensional and three- 
dimensional, have been examined and the exact solutions obtained by means of the 
Fourier, Mallin, and Hankel transforms. Until recently, the Fourier transform was 
the one most widely used. Filon's (ref. 7) investigation of the state of stress in an in- 
finite strip with a known exterior force was probably the first solution in the theory of 
elasticity solved by the use of Fourier integrals. Sneddon (ref. 8) investigated the 
elastic equilibrium of an infinite strip and, in particular, the case of a half plane. The 
application of the Fourier transform to the solution of the first fundamental problem for 
the half plane is discussed in the monograph of Novozhilov (ref , 9). With the use of the 
Fourier integrals in bipolar coordinates, the exact solutions of plane contact problems 
can be derived for a single punch and also for similar mixed problems of stresses in an 
infinite body weakened by a plane slit. The simplest of these problems was examined 
by Sneddon (ref. 10), and the solution was worked in Cartesian coordinates and led to a 
set of dual integral equations. A solution of the Boussinesq problem with the help of 
the Hankel transform is given in the books by Sneddon (ref. 8) and Tranter (ref. 11). 
Uflyand (ref. 12) obtained an exact solution of the Boussinesq problem by the use of the 
Papkovich-Neuber function of the mixed problem for an elastic layer, on one boundary 
of which are given the external forces and, on the other, the displacements. 
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Padfield and Sida (ref. 13) determined the elastic stresses that exist in a layer 
depressed symmetrically by two flat cylindrical punches. Muki (ref. 14) solved 
numerous asymmetric problems for a semi- infinite solid and a thick slab. Nelson 
(ref. 15) considered the problem of a thick plate loaded axisymmetrically. Greenwood 
(ref. 16) computed the stresses produced in the midplane of a slab by pressures applied 
symmetrically to its surfaces. The solution for the stresses in a plate that has been 
pressed between two spheres was given by Tu and Gazis (ref. 17). 

Many punch problems have been solved for the case when the punch is rigidly 
connected to the elastic body, when friction is present, or when two elastic bodies are 
in contact; but because they are not illustrated in this paper, they will not be mentioned 
here. The material cited in this section is not intended to be a complete bibliography 
on contact problems but rather a limited list of sources for information on specific con- 
tact problems. 


ELEMENTARY CONTACT PROBLEMS 


This section on elementary contact problems is presented as a natural starting 
point in the presentation of certain aspects of punch problems arising in the strength of 
materials. All simplifying assumptions of the theory of beams are taken into account. 
However, not all contact problems can be solved by this approach, and only the simpler 
cases are considered. The concept of contact plays an important role in the solution of 
problems involving the elementary theory of beams and plates because the concept 
changes the constraints of the problem. The problems presented here are the nonlinear 
bending of a cantilever beam, the nonlinear bending of a simply supported beam, and 
the problem of a thin plate resting on an elastic foundation. The solution to the last 
problem will be compared to the problem of a thick layer resting on an elastic founda- 
tion. 


The governing equation for the bending of beams is given by 


M_._-.l_ 

El " R 



( 1 ) 


or, because dy/dx is usually small, equation (1) can be approximated as 



Nonlinear Bending of a Cantilever Beam 


The problem is to determine the free-end (point B) deflection of a cantilever 
beam that is gradually forced into contact with a rigid circular surface (AC) (fig. 1). 
The solution to this problem may be found in reference 18 and is given as 


5 b- 6 i + Wfe 


(EX) 


( 3 ) 


where 6^ = the deflection of segment DB considered as a simple cantilever beam 

6 g = the deflection at point B caused by the slope at point D 

5„ = the deflection at point B representing the vertical distance of point D from 
the horizontal tangent at point A 


Nonlinear Bending of a Simply Supported Beam 

The problem is to determine the angle a , which defines the positions of the 
points of contact between the elastic beam and the rigid circular surfaces (fig. 2). For 
small angles, the result was presented in reference 18 as 


a = 


P g (L - 2Ra) 2 
16EI 


(4) 


Solutions for the problems involving the nonlinear bending of both the cantilever beam 
and the simply supported beam are obtained from the condition that at the points of 
contact the deflection curve is tangent to the supporting surfaces. 


Thin Plate Resting on an Elastic Foundation 

The problem is to determine the vertical displacement equation of a thin, infi- 
nitely extended plate that is resting on an elastic foundation and whose surface is sub- 
jected to a concentrated load (fig. 3). The reaction of the elastic foundation is 
assumed to be proportional to the deflection in the z-direction w and is given by Kw. 
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The constant K, expressed in units of pressure per unit deflection, is called the mod- 
ulus of the foundation. The general equation to be solved is in the form 



1 dVd 2 w 
* dr Adr 2 



q - Kw 
D 


(5) 


where q is the normal load distribution on the surface, h is half the plate thickness, 

and D is the plate rigidity given by D = 4Eh 2 /l2(l - y 2 ). For q equal to zero every- 
where except at the origin and also for the deflection vanishing for large values of r, 
the solution for equation (5) is given by 


-P 1 

w = "W kei(r/1) 


where 1 = (D/K)*^, and the subgrade reaction is 


-P Q kei(r/l) 


°z ~ .2 
2tt1 


In equation (7), kei(x) is defined by the relations 


Kq(xvCi) = ker(x) ± kei(x) 


( 6 ) 


(7) 


( 8 ) 


where Kq is a Bessel function of the second kind. This type of problem was originally 
solved by Hertz (ref. 19), and his data are presented in graphic form in figure 4. 


TWO-DIMENSIONAL PROBLEMS OF A HALF PLANE 


This section involves solutions for both the first fundamental problem and the 
punch problem in a two-dimensional half plane. However, before proceeding, a dis- 
tinction will be made between the first fundamental problem and the punch problem. 
The first fundamental problem involves the determination of the elastic equilibrium of 
a body when the pressure on its surface is a prescribed function (fig. 5). The second 
fundamental problem involves the determination of the elastic equilibrium of a body 
when the displacements on its surface are prescribed. The mixed boundary problem 
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is one in which the surface displacements are prescribed on only a portion of the bound- 
ary while the surface pressures are prescribed elsewhere. The problem of the punch 
is a special case of the mixed boundary problem in that the surface displacements are 
prescribed on a portion of the boundary, but the normal stresses are equal to zero ev- 
erywhere else (fig. 6). 

Two different but equivalent methods are used to solve the first fundamental prob- 
lem and the punch problem. The solutions derived by both the Cauchy integral and the 
integral transform will be listed preceding the appropriate section. The examples that 
are presented consist of a point force on the boundary, a uniformly distributed load on 
part of the boundary, a straight horizontal-base punch problem, and a wedge-base punch 
problem. 

The basic equations of the classical (infinitesimal) theory of elasticity that govern 
the problems solved here are given as follows in Cartesian coordinates. 

1. Equilibrium (with body forces neglected) 


do 


0CT 


_2L + _JEI = 0 

9x + 3y 


(9a) 


do 


xy 


do 


y _ 


J ^ s o 

3x + dy 


(9b) 


2. Compatibility 




.2 

9 y. 


xy 


3x 3y 


( 10 ) 


where e 


= 

x 3x ’ 


e = 

y 


3v 

ay’ 


and y. 


xy 


l/3u 

2\3y 



3. Stress- strain 


a 

x 



+ e y) + 2tie x 


(11a) 




(lib) 
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o = X(e +€\ + 2ae 
y V x yj Ty 


(He) 


First Fundamental Problem for Half Plane 

Let P(1 q) be the pressure distribution on the surface, and let denote the 

x-coordinate when y = 0 (fig. 7). The following formulas (ref. 10) are obtained by ap- 
plying Fourier transforms to the partial differential equations of the problem, by ob- 
serving boundary conditions, and by inverting the transform. 


cr 

X 




(x - t) 2 p(t)dt 



2 




(x - t)p(t)dt 



a = - 

y 



u = 



(x - t)y 

, + \2 2 
(x - t) + y 


dt 


v 


TTfl 



(1 - ^)(X - t) 

7 772 ? 

(x - t) + y 



where p 


,(t) = f p(t)dt. 
L J-OO 


( 12 ) 


(13) 


(14) 


(15a) 


(15b) 
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However, if the Cauchy integral method is used and the known equations (ref. 2) 
are utilized, the following formulas are obtained. 

°y * i<* xy = ®(z) - $(z) + (z - z )$’(z) (16) 

a x + dy = 2[Re®(z)] (17) 

2m |^ + i ^ = (3 - 4v)A(z) + $(z) - (z - z)l’(z) (18) 


where (’ ) denotes differentiation. 

For the first fundamental problem, the following equation is obtained. 


®(z) = 


X f p ( t o) ~ 'Wo 

2iri J L ( t o- z ) 


(19) 


where p(1q) and T^tq) are the normal and shear loading on the boundary, respec- 
tively. 

Concentrated force on the boundary . - For the case of a concentrated force acting 
on the boundary (fig. 8), the following formulas are obtained for the first fundamental 
problem of a half plane. 


a 

x 


2P(x - c) 2 y 




a 


y 



( 20 ) 


(21) 


( 22 ) 
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Equations (20), (21), and (22) can be obtained by the Cauchy integral equa- 
tions (16), (17), and (19) or by the integral transform equations (12), (13), (14), and 
(15). The solution to this problem is important because with it, the solution for any 
distribution of a load can be obtained by simple integration. 

Uniformly distributed load on a part of the boundary . - The solution for the prob- 
lem of a concentrated load on a boundary proves to be useful in solving the problem of 
the uniformly distributed load on a part of the boundary. The solution to the latter 
problem may be obtained by integration of equations (20), (21), and (22) with respect 
to x in the limits from -a to a, where -a and a are the extremes of the loading 
in the x-direction. The solution may also be obtained directly from equations (12) to 
(15) or from equations (16), (17), and (19). The solution is given as 


CT = 
X 






2pay(x 2 - y 2 - a 2 ) 


7T 


( 2 2 2 \ . 2 2 
[x + y - a ) + 4a y 






2 


a 2 ) 

, 2 2 
+ 4a y 


(23) 


(24) 




4paxy 2 



, 2 2 
+ 4a y 


(25) 


Punch Problems for the Half Plane 

In solving punch problems, as in solving the first fundamental problems, the 
basic equations (9), (10), and (11) must be satisfied, but a different approach is nec- 
essary because the boundary conditions are different. 

There are two classes of punch problems. The first class involves the type in 
which the surface of contact is known and does not change in the process of loading. 
The second is the class in which the surface of contact is to be found from the solution 
of the problem. A typical example of the first type is the flat-bottom punch (fig. 9). 

In the second type of punch problem, the surface contact depends on the proper- 
ties of the bodies in contact, on the profile of the punch, and on the applied force. The 
solution obtained by integral transform (ref. 8) is given by 


Vo) "Vo)/ Vo) 2 


(26) 


13 



where 


'( t o ) 1 


is the constant part of the deflection and 



( 27 ) 


The pressure under the punch is given as 



1 (n/2 - 1/2)’. a 

where A = — 7= s-r . By the Cauchy integral method (ref. 2 ), the solu- 

n ^ (l - * 2 )(n/2 - 1)'. 

tion for the second type of punch problem is also given by 


P(‘o) - 


M 


ttU - ^(to - a)(b - t 0 ) 


V(t - a)(b - t)f’(t) 

t - t„ 


dt + 


Wo ~ a )( b - t o) 


( 29 ) 


where a and b are the extreme points of contact, f’ (t) is the derivative of the punch 
profile, and Pq is the total vertical load. 

Cylinder indenting a semiplane. - The problem of a cylinder indenting a semiplane 
is an example of a punch problem in which the surface varies (fig. 10 ). The equation 
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for the profile of a cylinder is y = 2R = f^t^, and the contact pressure is, there- 
fore, found to be 


P 




(30) 


To find the value of a as a function of Pq, the contact pressure is observed to be zero 


at P(-a) and P(+a). Therefore, lim P(t n ) = 0, so that a = ± 

t_— a ' u/ 

0 

the pressure P^t^ is completely determined. 


P 0 R(k + l) 1 1/2 


2-nii 


. Now 


Wedge-shaped punch . - Another example of a punch with varying surface of con- 
tact is a wedge-shaped punch (fig. 11). The profile of the punch is given as 
f (to) = b + e^l - tg^/a) with 0 s tg < a. Here, the method of solution is exactly the 

same as for the last' problem, except that in this problem a corner exists at t^ = 0; 

therefore, the singularity appears at that point. The pressure under the punch is given 
as 


(31) 


In the classical (infinitesimal) theory of elasticity, the deflections and their par- 
tial derivatives are assumed to be small. Therefore, the squares and products of these 
quantities are neglected when deriving the basic elasticity equations . In this problem 
and in that of the flat-base punch, there are points on the contact surface where 
2/2 

9 v/ 0x (change in slope in the x-direction) is by no means negligible; therefore, at 
those points, singularities exist. 


FIRST FUNDAMENTAL AND PUNCH PROBLEMS 
FOR THE HALF SPACE 


The solutions to the first fundamental problem and the punch problem for a three- 
dimensional half space are presented in this section. The equations for these general 
solutions were derived by the use of the Hankel transform and are presented in the ap- 
propriate section. The first fundamental problems considered involve a concentrated 
force acting normal to the boundary and a uniform pressure distributed over a circle. 
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The punch problems that are presented involve a flat-base cylindrical punch and a 
conical-base punch. Cylindrical coordinates will be used throughout this section. 

The solution to the axisymmetrical stress problem must satisfy the following 
basic equations of classical theory of elasticity, as given in cylindrical coordinates. 

1. Equilibrium (with body forces neglected) 


2 . 


where y 
3. 


K 3lT rz „ 

dr + dz + r ' 0 

(32a) 

da da a 

rz + z + rz = 0 
dr + dz + r 

(32b) 

Compatibility 


2 

- fa - a e) + 1 l v far + + \) = 0 

(33a) 

v °e + r 2( a r ' °e) + 1 + v dr{ a r + (J 9 + a z) “ 0 

(33b) 

2d d d 2 

= ar 2 + rar + az 2 


Stress- strain 


„ f7s « \ 3 x] . 3w 

(34a) 


(34b) 

a z = x (fe + D u + (X + 2m) If 

: ; '• i ' -V: ‘ ' VV : : : • > V ' ' \ v : • V : i-V.- ■' • • ' O ,■ ' ;! , : : ' . : ' ^ ; 

(34c) 

/\ 3 X + 2 u\ . 3w 

a. = X -y + =•) u + A 

6 \ dr r J dz 

(34d) 
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The Hankel transform is used extensively in these sections. Therefore, it 
becomes necessary to define this transform and its inversion formula. The Hankel 
transform of the function f(r) given on the positive real line 0 < r < 00 is defined by 
the integral 


*„(!') -jf f(r)«y£r)r dr (35) 

where 0 £ £ < 00 , v > -1/2 and J^(x) is the Bessel function of the first kind. If the 
function f(r) is continuous in parts in any finite interval belonging to the region (0, 00 ) 
and if the integral J^°°|f(r)|\/r dr converges, then the Hankel transform exists. 

If the function f(r) also satisfies the Dirichlet conditions in any open interval 
0 < r < s, the Hankel inversion formula is 


f(r) = 



(4)Jj,(£r)£ d? 


(36) 


where 0 < r < 00 . 


First Fundamental Problem for a Half Space 

By first applying the Hankel transform to the basic elasticity equations, then ap- 
plying the appropriate boundary conditions, and finally using the inversion formula, the 
following solution for the first fundamental problem of a half space (ref. 10) can be ob- 
tained. The problem is illustrated in figure 12. 

If p(r) is the axisymmetric pressure distribution on the surface z = 0 and if 
the shear stresses are zero on this surface, then the equations for the solution of the 
problem in the half space are given by the following conditions and by equation (37). If 

f(r) = -p(r)/2ju and the total load P Q = ri(r)dr and Tq(^) = Hg[f(r); r — £] where 

H.( ) denotes the ith order Hankel transform, then 


(r, z) = -H^td - 2v) - |z] r 1 f(«)e" €z ;S - r} 

(37a) 

5 (r, z) = H 0 {[2(1 - v) + Sz].r 1 1(|)e“ 5z ;| - r| 

(37b) 
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CT rz (r,z) = -2 J LtzH 1 |jf (g)e _ ^ z ;4 - r] 


a z (r, z) = -2 J uH 0 [(l + £z)f (z)e~^ z ;£ - r] 


(37d) 


a r (r,z) = -2 M H 0 [Tl - £z)f (Oe"^ z ;£ hJu - 2 ^ _1 - zf(0e"5 ( 3 ?e) 

a 0 (r, z) = -4 4 ^H 0 [r(4)e" ?z ; * - ^ ~ ^ H 1 |jjX - 2v)f l - ^ f (|)e' |z ; ? - rj (37f) 

Concentrated force acting normal to a half space . - As an example of the first 
fundamental problem, the problem of a concentrated force acting normal to a half space 
(Boussinesq problem) will be presented. The problem is illustrated in figure 13. Be- 
cause the force is localized at the origin, the expression f(r) = k may be used 
where 6(r) is the Dirac delta function. 

Choose k so that the equation Pq = 4 rf(.r)dr is satisfied. Note that 

p 0 p 0 6(r) _ P Q 

k * 4i/7 and ’ therefore > f ( r ) = 47r/xr »' thus, f 0 (|) = 4 From the preceding equa- 
tions, it can be shown that 


*0 ) (1 - 2v) |7 z 

U r _ 477 > r 7T TJ72 

(r 2 + z 2 ) 


' 2 2 ' 
*.r + z , 


(38a) 


P 0 2(1 - v) , z 2 

W — -j — ' — — ■ . . — 4* — - — — — — 

Z 4 TT/x f o 2 / 2 

V r + z l (r + z 


(38b) 


a (r, z) = 
rz v ? 


3 p 0 z 2r 


27r(r 2 + z 2 ) 


a z (r,z) = 


3P 0 Z ' 


/ 2 2V 
r(r + z ) 
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z) = 


z 


(38e) 


P 0 (l - 2v) 

z 1 

2ir 

(2 2^ 2 r 2 

(r +z ) 




CT r ( r > z) 




3r 2 z 


572 ' 


(r 2 + z 2 ) 


(38f) 


Uniform pressure over a circle . - The solution for a uniform pressure distributed 
over a circle may be obtained by simple integration of the concentrated force solution 
(eq. (38)) or by the application of equation (37). This problem is illustrated in fig- 
ure 14. The case of the surface pressure distributed uniformly over a circle of radius 
a is represented by 


p(r) = -pH(a - r) 


(39) 


where H(co) = 1 when co > 0 and H(o>) = 0 when u> < 0. The total 


load 


to 7ra p, or p = 


f(r) = 


P Q H(a - r) 

27ra 2 /j. 


P 0 p 0 
— 2 ? and equation (39) becomes p(r) = 2 H ( a “ r )* 

ra nar 

PqJj^) 


and f Q (£) - 2jrajljl | 


Pq is equal 
Therefore, 


If the following notation is introduced I (if/, $, y) 


f J^aOJ^r^^e'^d?, then 


u r (r, z) = - 1 - 2^)I(1, 1, -1) - zl(l, 1, 0) 

(40a) 

w z (r, z) = 27rajU 2(1 - v)I(l, 0, -1) + zl(l, 0, 0) 

(40b) 

<I rz (r " z)= ‘I! 1 * 1 ' 0 ’ *> 

(40c) 

°zz “ ‘ 1 ; I(1> °- 0) + zl(1 > °> '> 

(40d) 
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v = ■ m I(1 > °> 0) ■ zl(1 > °’ 1} + ii < 1 - *> -*> - *> °) ( 40e ) 

I(1 > °> - JiF (1 ' 2v)1 ^’ -D - i» °) ( 40f ) 


Numerical values for I(»^, /3, y) may be found in reference 20. 


Punch Problem for the Half Space 

When solving the punch problem for the half space, as in solving the first funda- 
mental problem, the partial differential equations of elasticity (eqs. (32), (33), and 
(34)) must be satisfied. The boundary condition for the punch problem is a known dis- 
placement under the punch and zero normal stress elsewhere on the boundary, as illus- 
trated in figure 15, By transforming the equations, using the boundary condition, and 
inverting the results, Sneddon obtained the following solution for the axisymmetric 
punch problem (ref. 10). If the displacement under the punch is w(r, 0), then 


w(r, 0) = b + f(r/a) (41) 

where 0 < r < a, b is as yet unspecified, and f(r/a) is a prescribed f miction; also, 

a z (r,0) = 0 (42) 


where r > a. 

Define the functions y(t) and such that 



where x = r/a, and 


m = x(t)cos (£t)dt (44) 

•'o 
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Then the total Pq is given as 


P 0 - a'(i ? ivi / L rJ 0 (?r/a)dr 


(45) 


The pressure directly under the punch is given as 


% <r ’ 0) “ afT^T) 


*( 1 ). 


vr 


-jfe 


: dt 


X 2 "X 


(46) 


When the pressure under the punch is known, the problem may be solved by the method 
given in the first part of this section. 

Flat-base cylindrical punch . - As an example Of the type of punch whose surface 

area does not vary with load application, the punch in the form of a flat-base cylinder 

(fig. 16) will be considered. For the horizontal-base punch, w (r, 0) = e where 

z 

0 ^ r < a. By equations (41) and (42), b = e and f(r/a) = 0. Therefore, by equa- 

2e 

tion (43), it can be seen that x(t) = — ; and by equations (44) and (45), the total load on 


the punch is Pq - ^ . From this, the equality e = — is obtained. The 

pressure under the punch is given by equation (46) as 


p 0 d - -> 


-' 2 (r. 0) = 


2jue 


ir(l 


- v)Va 2 - r 2 


(47) 


If the expression for 6 is substituted into this last equation, the pressure under the 
punch will be obtained as a function of the applied load. 


tf z (r, 0) = 



(48) 


Conical-base punch . - An example of the punch problem in which the surface of 
contact varies is the conical-base punch, illustrated in figure 17. For this problem, 

w (r, 0) = b + e - ^ where r < a. Let m = r/a; then w (r, 0) = b + e (1 - m), and 

Z 2i ■% . Z 
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therefore f(m) = e (1 - m). The total load Pq as a function of e then becomes 


and 



4 u are 


b = -f(0) - |L(s)*E 

Jy/TJ 

b = e(|-l) 

*<*>■* + f( 1 -f t ) 


(49) 

(50) 

(51) 

(52) 


X’(t) = -e 


a z ( r, 0) = - cosh 1 (a/r) 


where 0 < r < a; or equivalently 


CT z (r, 0) = - — 5.- cosh" 1 (a/r ) 
4a 7r 


(53) 

(54) 


(55) 


where 0 < r < a. Note the similarity between the curves for the pressure under the 
punches of the two-dimensional case and those of the three-dimensional axisymmetric 
case. 


FIRST FUNDAMENTAL AND PUNCH PROBLEMS FOR AN ELASTIC 
LAYER RESTING ON A RIGID FOUNDATION 


In this section, the solutions are sought for the first fundamental and the punch 
problems for an elastic layer resting on a rigid foundation. The basic configuration 
for this problem is shown in figure 18. The equations of the classical theory of elas- 
ticity (eqs. (32), (33), and (34)) must be satisfied by the solution for an elastic layer. 
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The solutions presented in this section were derived by the use of Hankel transforms. 
Two examples are presented for the elastic layer resting on a rigid foundation. The 
first example is that of a uniform pressure distributed over a circle, and the second 
example consists of a flat-ended cylindrical punch being forced into the free surface of 
the layer. 


First Fundamental Problem for an Elastic Layer 
Resting on a Rigid Foundation 

For the first fundamental problem, the Hankel transform is applied to the basic 
elasticity equation, the appropriate boundary conditions are used, and the resulting ex- 
pressions are inverted to give the following solution for the elastic layer (ref. 21). If 
an elastic slab is resting on a rigid foundation with the surface loading axisymmetric 
about z equal to p(r), then the following equations may be used to solve the problem. 

Let m = r/b, p(r) = -2/if(r/b), cr^ = 0, and let Fq( w) be the notation for the 
zeroth order Hankel transform of f(r/b). Then 


2 — 

sinh wf n (w)J n (mw) 

w (mb, b) = -4(1 - u)b I — =— ; /» v — dw 

z' ’ ' ' ' Iq 2w + sin (2w) 


(56a) 


cr z (mb, 0) = -4/1 


[si 


[sinh(w) + w cosh(w)] J n (mw) 


2w + sinh(2w) 


dw 


(56b) 


(°r + a B + c z) 


z=0 


J r °°wfQ(w)sinh wJg(mw) 
n 2w + sinh(2w) 


(56c) 




f 

J 0 

J r OO 

n 


wf 0 (w)[(l - 2v)sinh(w) - w cosh(w)] JQ(mw) 


2w + sinh(2w) 


dw 


(w)[(l - 2v)sinh(w) - w cosh(w)] J^(mw) 
2w + sinh(2w) 


dw (56d) 


For the case of a load uniformly distributed over a circle (fig. 19), if H(a>) = 1 
when u> > 0 and H(u>) = 0 when w < 0 (where a> is any argument of H), then 
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p(r) = -p/jra^H(a - r) and 


p(r) = cr (r, b) = - PH (%~ — = -2/iif(r/a) 


77a 


pma/h-m^ _ PtJJtw) 

Let m = r/b; then f(m) = — - ■■■■ ~ so that f n (w) = s where r = a/b. 

27ra Z /j, 27ra pw 

stituting Fq(w) into the list of equations, 


<r r =J ^P'l- G l- F 2 + (1 - 2 ‘' )G 2l 

jU,77a ' — — 1 


a 0 = 


2Pt 


lita. 


- F 2 * <> ' 


a = - 
z 


2Pt(F 1 + .GjJ 

2 

juira 


where F^, G^, Fg, Gg denote the integrals defined by the equations 


'1 (m ’ T)= /„ 


w cosh(w)jQ(mw)Jj(rw) 
2w + sinh(2w) 


dw 


f sinh wJ 0 (rr 

¥ m > T)= j„ 2wTii 


(mw)J 1 (rw) 

sinh(2w) 


dw 


^ r cosh wJ 1 (mw)Jj(Tw) 

t) - — J 2w + sinh(2w) dw 


G 


/* OO 

/ sinh wJ ^(mw)J^(rw) 
j(m, T ) - ^ J 2w + sinh(2w)(w) dw 


(57) 

Sub- 

(58a) 

(58b) 

(58c) 

(59a) 

(59b) 

(59c) 

(59d) 
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Numerical values for equation (59) for various values of m and r are given in 
reference 22. Data for Sneddon’ s solution for contact pressure between a layer and the 
rigid foundation are presented in figure 20. 

Punch Problem for an Elastic Layer Resting 
on a Rigid Foundation 

When 0 ^ r < a, the equation u (r,b) = -D - w n (r/a) applies to the flat-ended 

z u 

cylindrical punch as shown in figure 21. The depth of penetration is D = e, Wq(x) - 1. 
It can then be shown that x (t) = ^ g(t) where g(t) represents the solution of the 

Fredholm equation. 


g(t) - / g(u)K(t, u)du 



The integral equation (60) has been solved numerically by Lebedev and Uflyand 
(ref. 23). The ratio 
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has been computed for a range of values of the ratio a/b; the results are given in the 
following table. 


a/b 

0 

0,5 

1.0 

1.5 

2.0 

P 0 (l - V) 
4a/ie 

1.0 

1.51 

2.20 

2.95 

3. 72 


SOLUTION FOR THE PROBLEM OF A THICK LAYER RESTING 
ON AN ELASTIC (WINKLER TYPE) FOUNDATION 


This section concerns the determination of the stresses and displacements in a 
thick layer that is resting on an elastic (Winkler type) foundation and whose surface is 
loaded by any known axisymmetric pressure distribution (fig. 22). This problem is of 
practical importance in the design of airport runways, building foundations, seals, con- 
crete roads, or any arbitrarily thick layer resting on an elastic foundation. Because of 
the generality of this problem, its solution also may be applied to the design of layers 
resting on rigid foundations and to the design of very thick layers that approach the half 
space. 

The method used in the solution of this problem is that of the Hankel transform 
which consists of transforming both the boundary conditions and the governing partial 
differential equations. The transformation of the partial differential equations reduces 
them to ordinary differential equations which, along with the transformed boundary con- 
ditions, may be solved by elementary means. Thus, the expressions for the transforms 
of the desired solution can be obtained, and by inversion of the transforms, the solution 
itself can be obtained. 

Considering bending only, Hertz (ref. 19) solved the problem of a thin, infinitely 
extended plate resting on an elastic foundation. Sneddon (ref. 21) obtained the solution 
to the problem of a thick layer while considering all the components of stress, but he 
considered the layer to be resting on a rigid foundation. It is believed that the problem 
solved in this section is of a more general nature than either of the previously men- 
tioned problems and, therefore, will reduce to Hertz's solution for small plate thick- 
nesses and to Sneddon’s solution for an infinitely stiff foundation modulus, 

The Winkler type of foundation is one in which the contact pressure, which exists 

between a rigid punch and the foundation, is proportional to the surface displacement of 

the foundation a = Kw. No confusion should arise between the Winkler foundation and 
z 

the semispace foundation because the contact pressure for the latter is given in the 
complicated expression of equation (46). 

The problem solution of a layer resting on a Winkler type foundation must satisfy 
the equations of classical elasticity (eqs. (32), (33), and (34)), For the case of sym- 
metry, the displacement components are (u, 0, w), and the nonvanishing stress com- 
ponents are cr^, o^, a , and o . The boundary conditions for this problem are 
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cr = 0 and cr = -p(r), where z = +h; and a =0 and cr = Kw, where z = -h. The 
rz z ” ’ rz z V 

Navier equations reduce to the following forms. 


2(1 - u)\ 



1 9u 
+ r 9r 



+ (1 - 2u) 


a 2 u a 2 w 
+ drdz 


= 0 


(63) 


(1 - 2v) 




+ 2(1 - v) 


a 2 w 

^7 



(64) 


The first order Hankel transform of the radial component of displacement will be de- 
fined as 


U 1<5’ Z > = I 


ru(r, z)Jj(£r)dr 


(65) 


Similarly, the zeroth order Hankel transform of the z component of displacement will 
be defined as 



rw(r,z)J Q (£r)dr 


(66) 


By first multiplying both sides of equation (63) by rJ^(£r) and then integrating 

with respect to r from zero to infinity, and finally using the well-known properties 
(ref. 8) of Hankel transforms, the original equation becomes 


u - 


2i,)D - 2(1 




£Dw q = o 


(67) 


where D = Also, by multiplying equation (64) by rJ^r) and then integrating 
with respect to r from zero to infinity, the original equation becomes 



v)D - (1 - 2u)£ 


Wq + ^Du 1 


0 


( 68 ) 
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The solution to equations (67) and (68) will be taken in the following form. 

u x = sinh(£z) + A 2 cosh(£z) + AgZ£ sinh(^z) + A^z£ cosh(^z) (69) 

Wq = Bj sinh(£z) + B 2 cosh(£z) + B^zi; sinh(|z) + B^zt, cosh(^z) (70) 

After substituting these expressions into equations (67) and (68) and noting that for any 
z they must be true, the following relation is found among the constants: 


CO 

m 

i 

n 

<j 

(71a) 

A 3 " " B 4 

(71b) 

A 2 = -Bj + (-3 -f 4y)B 4 

(71c) 

A^ = -B 2 + (-3 + 4i')B 3 

(71d) 


Therefore, the transform of u and w may be expressed in terms of only the B’ s as 


u^ = -B 2 + (-3 + 4 i^)B 3 sinh(4z) + -Bj + (-3 + 4y)B 4 cosh(£z) 

- B 4 z£ sinh(^z) - BgZ^ cosh(ijz) (72) 

w Q = Bj sinh(|z) + B 2 cosh(£z) + B 3 z£ sinh(^z) + B^zE, cosh(£z) (73) 

The shear stress expressed in terms of displacement is + |^. 

The first order Hankel transform of o yields 

rz 

V = I r<T rz J l ( ^ r)dr = m(d5i - ew 0 ) (74) 
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Therefore, 


°r Z = M^{-[T3 - 4i^)B 3 + B^jcosh(|z) - [73 - 4i^)B 4 + B^]sinh(£z) 

- B 4 [sinh(?z) + z£ cosh(£z)] - Bg[cosh(^z) + z£ sinh(£z)] 

- B^ sinh(^z) - Bg cosh(^z) - B^z^ sinh(^z) - B^z^ cosh(£z)j (75) 


The boundary condition on o at z = ± h is that cr =0 and, therefore, a =0 
J rz rz ’ rz 

also. The following equations are obtained by using the latter boundary condition and 

equation (75). 


[4(1 - ^)Bg + 2B^j cosh(£h) + JT(1 - y)B 4 + 2B^J sinh(£h) 
+ 2B 4 h£ cosh(£h) + 2B 3 £h sinh(£h) = 0 

[4(1 - y)Bg + 2B^jcosh(£h) - [4(1 - v)B^ + 2B^j sinh(£h) 
- 2B 4 h£ cosh(£h) + 2B 3 £h sinh(£h) = 0 


(76) 


(77) 


Addition of equations (76) and (77) yields 


[2(1 - y)B 3 + B^jcosh(£h) + Bg^h sinh(£h) = 0 


Subtracting equation (77) from equation (76) results in the following equation. 


j^2(l - y)B 4 + B^j sinh(£h) + B^h cosh(ijh) = 0 


(78) 


(79) 


From equations (78) and (79), the following relationships among the B’s can be ob- 
tained. 


B 2 = -B 3 [2(l -v) + £h tanh(£h)] 
Bj = -B 4 [2(l - v) + £h coth((jh)] 


(80) 

(81) 
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as 


The expressions for Uj and w Q may be written in terms of Bg and B 4 
follows. 


Uj = [2v - 1 + £h tanh(^h)]Bg sinh(|z) + [2v - 1 + |h coth(£h)]B 4 cosh(£z) 

-B^z sinh(^z) - Bgi;z cosh(^z) (82) 

w 0 = -[2(1 - u) + £h coth(|h)] B 4 sinh(£z) - [2(1 - u) + |h tanh(^h)] Bg cosh(£z) 

+ Bg^z sinh(i-z) + B^z cosh(iz) (83) 

Also, 


Therefore, 


a - 2 mU - ^)dw 
z (1 - 2vjdz 


+ 


_^jiu_fdu l\ 
1 - 2i'\3r + rj 


u 



2ju 


(1 


^)w 0 ’ + 


1 - 2u 



(84) 


(85) 


When u 


j and Wq are substituted into equation (85), a z becomes 

* = -2 m ||[1 + £h coth(^h)]B 4 cosh(|z) + [1 + £h tanh(£h)]Bg sinh(£z) 

-4BgZ cosh(ijz) - £B 4 z sinh(^z)^ 


( 86 ) 


The boundary conditions on a are as follows: where z = +h, o = -p(r) and, there- 

z z 

fore, cr = -P * and where z - -h, cr = Kw and, therefore, cr = Kw n . Here, P n 

2 g w Z Zq U U 

is the zeroth order Hankel transform of p(r). By applying the last expression for a 

Z 0 
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(eq. (86)) to the boundary conditions, the following two equations are obtained. 

P Q = 2/iijj[l + £h coth(£h)] cosh(£h)B 4 + [1 + £h tanh(£h)] sinh(£h)B 3 

- £h cosh(^h)B 3 - £h sinh(£h)B 4 j (87) 

k|[2(1 - v) + £h coth(£h)]sinh(£h)B 4 - [2(1 - v) + |h tanh(£h)] 
x cosh(^h)B 3 + ijh sinh(£h)B 3 - £h cosh(£h)B 4 j = 

- -2jtx?{[l + £h coth(£h)] B 4 - [1 + |h tanh(^h)] sinh(£h)B 3 

+ £h cosh(^h)B 3 - £h sinh(£h)B 4 j (88) 


Equations (87) and (88) are two equations which can be solved simultaneously so 
that the constants B 3 and B 4 may be determined as 


_ ~ P 0 I [2K(1 - y)sinh 2 (gh) + M g[sinh(2gh) + 2gh] [cosh(gh) 
2 ^ | 2ju^ _(2i;h) 2 - sinh 2 (2£h7] - K(1 - y)[4£h + sinh(4£h)] 


(89) 




{2K(1 - y)cosh 2 (gh) + jug[sinh(2gh) - 2^h] } sinh(gh) 
2jU^jj2^h) 2 - sinh 2 (2£hj] - K(1 - i/)[4£h + sinh(4£h)] 


(90) 


From the previous expressions for the relation among the B’s, values for Bg and B 3 
are determined as 


2K(1 - y)sinh 2 (gh) + jn4[sinh(2gh) + 2gh] 

2/x^ [^(2^h) 2 - sinh 2 (2£h)] - K(1 - f)[4£h + sinh(4£h)] 

x [2(1 - y)cosh(^h) + £h sinh(£h)]| 


B -5° 
2 ~ M £ 


(91) 
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B = ^2 

B 1 & 


2K(1 - i/)cosn (gh) + M £[sinh(2£h) - 2|h] 


2p£[(2£h) 2 - sinh 2 (2ijh)] - K(1 - v)[^h + sinh(4£h)] 


x [2(1 - y)sinh(£h) + £h cosh(£h)J 


( 92 ) 


Recalling the relation among the A’s and the B's, the following expressions are 
obtained for the A' s . 


- p o! 

|2K(1 - y)sinh 2 (£h) + ju£[sinh(2£h) + 2£h] } cosh(£h) 

| 

2<it[(2£h) 2 - sinh 2 (2£hl] - K(1 - ^)[4^h + sinh(4£h)] 


x [2v - 1 + |h tanh(^h)] 


(93) 


- p ol 

|2K(1 - y)cosh 2 (£h) + [sinh(2£h) - 2£h] [sinh(^h) 

j 

2 M ^[(24h) 2 - sinh 2 (2£h7] - K(1 - y)[4£h + sinh(4£h)] 


x [2u - 1 + £h coth(£h)] 


(94) 


A -!• 

3” Ml 


2K(1 - v ) cosh (^h) + M £[sinh(2gh) - 2£h] 


2jLt^jT2^h) 2 - sinh 2 (2£h)] - K(1 - v)[4:%h + sinh(4^h)] 


sinh(£h) (95) 


A 3 
4 Ml 


2K(1 - y)sinlT(gh) + ^J[sinh(2|h) + 2gh] 


2 jU ^[T24h) 2 - sinh 2 (2£h)J - K(1 - i/)[45h + sinh(4£h)] 


cosh(£h) (96) 
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By substituting equations (89) to (96) into equations (82), (83), (84), and (86), the fol- 
lowing expressions are obtained for u 1 , w n , a , and a 

i u Zq rzj 


Uj = - {[2v - 1 + ijh tanh(i-h)] sinh(£z) - £z cosh(£z)} 
P 


x 


Ml 


{2K(1 - y)sinh^(gh) + ^g|sinh(2^h) + 2|h] }cosh(gh) 
2/Lif[(2|h) 2 - sinh 2 (2|h)] - K(1 - v)[4ijh + sinh(4£h)] 


-\[2v - 1 + £h coth(|h)]cosh(iz) - i-z sinh(£z)| 


m? 


|2K(1 - i/)cosh Z (gh) + Mf[sinh(2gh) - 2gh{f sinh(£h) 


2/i« 


(2£h) 2 - sinh 2 


- K(1 - i/)[4£h + sinh(4i-h)] 


(97) 


w n = -{-[ 2(1 - u) + |h coth(ijh)] sinh(i-z) + i-z eosh(|z)} 


' mT 


{2K(1 - y)co.sh (gh) + juijfsinh(2ijh) - 2i;h] } sinh(^h) 


2 ni 


(2£h) 2 - sinh 2 


K(1 - *»)[4$h + sinh(4ijh)] 


-{-[2(1 - u) + ijh tanh(i-h)]cosh(ijz) + i-z sinh(i-z)} 


Ml 


{2K(1 - v)ginh 2 (£h) + u g[slrih(2gh) + 2gh] [ cosh(gh) 


2^1 [(2|h)^ - sinh 2 (2£h^ - K(1 - y)[4i;h + sinh(4i;h)] 


(98) 


a = 2{[1 + ijh coth(i-h)] cosh(^z) - i-z sinh(i-z)} 

Z 0 

P Q |2K(1 - v)cosh 2 (?h) + /ig[sinh(2gh) - 2£h]|sinh(£h) 
2jn||T2|h) 2 - sinh 2 (2£h)] - K(1 - i/)[4gh + sinh(4£h)] 
+ {[1 + £h tanh(i-h)]sinh(i-z) - £z cosh(ijz)[ 

P Q |2K(1 - v)sin» 2 (|h) + /ii[sinh(2£h) + 2£h]}cosh(ih) 
2 j ug[(2£h) 2 - sinh 2 (24hJl - K(1 - v)[4^h + sinh(4£h)] 


(99) 
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a = -|2£[h cosh(£z)tanh(£h) - z sinh(^z)] } 

1 

P Q |2K(1 - i/)sinh 2 (4h) + ju£[sinh(2£h) + 2|h]Jeosh(£h) 
2/i^(T 2£h) 2 - sinh 2 (2^h7J - K(1 - v)[ 4£h + sinh(4£h)J 
+ |2£[h coth(£h)sinh(£z) - z cosh(£z)]} 


x 


P q | 2K(1 - t')cosh 2 (|h) + ju£[sinh(2£h) - 2£h]|sinh(£h) 
2/ig [(2^h) 2 - sinh 2 (2£h7| - K(1 - M)[4^h + sinh(4|h)] 


(100) 


The preceding expressions for the Hankel transforms of u, w, a , and o 

z rz 

only need be inverted to obtain the desired solutions. The Hankel inversion formula is 



ij,t, z)J^r)S d? 


( 101 ) 


where 0 < r < 00 . Therefore, 


u(r, z) = j UjJj^r)? d£ 


( 102 ) 


w(r, z) = J w Q J 0 (£r)£ d| 


(103) 


or 

z 



0 (|r)S d£ 


(104) 


a 

rz 



J r (?r)? d£ 

1 1 


(105) 


Substituting equations (97) to (100) into the appropriate inversion formula yields the 
desired solution in terms of any prescribed surface pressure p(r). 
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u(r, z) - - f \[2v - 1 + £h tanh(4h)j sinh(£z) - |z cosh(£z)f 

J r\ 


Jj(€r)P 0 . 


2K(1 - v)sinh 2 (£h) + /i£[sinh(2£h) + 2£h] jcosh(£h)d£ 


m{ 2 M |[(2£h) 2 - sinh 2 (24hT| - K(1 - i>)[4£h + sinh(4£h)]}- 

r 

I {[2^ - 1 + £h coth(£h)]cosh(£z) - £z sinh(£z)[ 

•Vi 


X 


J^iOP^KU - i/)cosh 2 (?h) + M £[sinh(2£h) - 2£h]}sinh(£h)d£ 
/u.j2^£ [(2£h) 2 - sinh 2 (2£h)] - K(1 - y)[4£h + sinh(4£h)]| 


(106) 


w(r, z ) = - I 

Vrt 


{-[2(1 - v) + £h coth(£h)]sinh(£z) + §z cosh(£z)[ 


J 0 (^r)P 0 «j2K(l - v)eosh 2 (£h) + M |[sinh(2£h) - 2£h]}sinh(£h)d£ 
M |2 M ^|72^h) 2 - sinh 2 (2£h)] - K(1 - v)[^b. + sinh(4£h)]| 

f {‘-[2(1 - v) + £h tanh(£h)]cosh(£z) + £z sinh(|z)[ 

0 


J Q (^r)P 0 |2K(l - ^)sinh 2 (|h) + [sinh(2£h) + 2£h] jcosh(£h)d£ 
ji^f^h) 2 - sinh 2 (2£h]] - K(1 - u)[A^h + sinh(44h)]} 


(107) 


a z (r, z) = 2J | [1 + £h coth(£h)]cosh(£z) - £z sinh(£z)f 


J 0 (£r)P 0 £j2K(l - i/)eosh 2 (£h) + M i[sinh(2£h) - 2£h]Jsinh(£h)d£ 


+ 2 J |[1 .+'ih tanh(£h)]sinh(£z) - £z cosh(|z)} 


J 0 (?r)P 0 ^2K(l - y)sinh 2 (£h) + M |[sinh(2|h) + 2£h]lcosh(£h)d£ 

X p— " n — n — j “““ (108) 

2 M 4[(_24h) - sinh (2£h]J - K(1 - i/)[4£h + sinh(4£h)] 
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ov (r, z) = - f [2i;h cosh(£z)tanh(£h) - 2£z sinh(£z)j 
rz j q 

J^r^P^Kd - „)sinh 2 (£h) + ^|[sinh(2|h) + 2|h]}cosh($h)d$ 
2/i? [(2£h) 2 - sinh 2 (2^hj| - K(1 - i^h + sinh(4£h)] 

- f [2£h coth(£h)sinh(i;z) - 2£z cosh(ijz)] 

J 0 

JjCir^PQ^Kd - y)cosh 2 dh) + M 4[sinh(2|h) - 2£h])>sinhdh)d£ 
2^ij [?2^h) 2 - sinh 2 (2£hTJ - K(1 - v)[4£h + sinh(4£h)] 


(109) 


NUMERICAL RESULTS ON THE CALCULATION OF CONTACT STRESS 
BETWEEN A LAYER AND ITS ELASTIC FOUNDATION 


In the preceding section, equations were derived for displacements and stresses 
in an elastic layer resting on an elastic foundation and subject to any prescribed axi- 
symmetric surface pressure. A comparison will be made between these equations and 
those derived by Sneddon (ref. 21) for a layer resting on a rigid foundation. The equa- 
tions derived in the last section will be used to determine the distribution of contact 
pressure between a layer and its foundation for the case of the surface pressure uni- 
formly distributed over a circle. The equations derived in the last section involve in- 
finite integrals which, for the purpose of plotting the results, were numerically 
integrated. Numerous examples are solved for different foundation moduli and differ- 
ent circle radii. The results are illustrated in figure 23. 

The following equation was presented by Sneddon (ref. 21) for the contact stress 
between an elastic layer and its rigid foundation. 


CT 

Z 




[tj cosh(T}) + sinh(r?)]T 7 P(r//d)J 0 (pr 7 )dT 7 
~ 2r\ + sinh(2r?) — — 


(HO) 


where d is the layer thickness and p = r/d. The equation for the stress in the 
z-direction for a thick layer resting on an elastic foundation was solved in the preceding 
section and was given in equation (108). The following expression can be obtained for 
lim cr (r, z) by dividing both the numerator and the denominator by K(1 - u) and taking 
K— oo z 
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the limit as K approaches infinity 


lim or (r, z) = -4 f J J£r)P^|ili-£ 
K-°° Z •'o u u ( 


l;h coth(£h)]cosh(i;z) - gz sinh(gz)} 
4^ h + sinh(4|h) 


x cosn (£ h)sinh(£ h) 


|[1 + |h tanh(|h)]sinh(£z) - £z cosh(£z)}sinh^(£h)cosh(£h)) At 
+ 4|h + slnh(^h) I** 5 


( 111 ) 


The equation for contact stress for z = -h is given by the preceding equation. 


lim'cr (r, -h) 

K-00 z 


f^T /r r yp ,fsinh(2^h) + (2|h)cosh(2?hT| , fc 
| q J Q ^r)P 0 ^ 4 ^h + sinh(4£h) J d| 


( 112 ) 


Let rj = 2£h; therefore, £ = 77/2I1 and d£ = dr]/2h. Also, let p = r/2h and = a/2h. 
Then 


lim a (r, -h) = 
K-°° z 


lim a (r, -h) = - 

K-co z 


, p o” 

sinh(r]) + 77 cosh (r]) 

/ 2h 

_ 2 ?] + sinh(2rj) 


1 ft ^ r \p n sinh ( r J) +7 ? cosh(7j) , 
~2J Q J OP P 0^[_“2}j + sinh(2n)-^J a7? 


(113) 

(114) 


Because h is only half the layer thickness [h = (l/2)d] , the following equation 
is obtained in terms of d. 


lim o (r, -h) 

K-00 z 


f 


d^O 


J 0 (P*?)P 


.(7?/d)7? 


s 


sinh(r?) + r] cosh ( 77 ) 
2 tj + sinh(2rj) 




dr? 


(115) 


This solution is the same as the one obtained by Sneddon. 
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Derivation of the Contact Equation 

In this section, curves for the contact pressure between the layer and its elastic 
foundation will be given for a surface loading of a uniform pressure distributed over a 
circle (fig. 23). The following material presents a step-by-step derivation of the con- 
tact pressure equation to be numerically integrated. 

By using the expression for a (r, z) in equation (108) and evaluating it for 

z 

z = -h, the following equation is obtained. 


CT z (r, -h) = 2 


r 


J 0 ttr) p 0 5 


[sinh(^h)cosh(^h) + gh] {2K(l - u)cosh 2 (gh) + jug[sinh(2gh) - 2gh][ 


2/i|[(2gh) 2 - sinh 2 (2£h)] - K(1 - i/)[4$h + sinh(4£h)] 
[-sinh(£h)cosh(|h) + £h]| 2K(1 - v)sinh 2 (£h) + pi[sinh(2(;h) + 2|h]} 

"J" ' ‘ f— ■“ K n — i ‘ “ 1 — — 


d^ 


2ju^ |j24h) 2 - sinh 2 (2£h7] * K(1 - y)[4^h + sinh(4|h)] 


(116) 


When terms are collected and the hyperbolic identities are used, equation (116) reduces 
to 


<* z (r, -h) 



J 0 (£r)P 0 £ {K(l - v)[sinh(2£h) + 2£h cosh(2£h)]}d£ 
p£[(2(;h) 2 - sinh (2£h)] - K(1 - i^)[2^h + sinh(2|h)cosh(2£h)] 


(117) 


If p(r) is described as p(r) = p where 0 < r < a and p(r) = 0 where r > a, 

then P q = j Jj(a£)p. By inserting this into equation (117), the following expression 

for ct (r, -h) is obtained, 
z 


a z (r, -h) = ap 



J 0 (|r)J 1 (a^){K(l - ^)[sinh(2|h) + 2fh cosh(2£h)]f d£ 
p|[( 2£h) 2 - sinh 2 (2£h)] - K(1 - ^)[2^h + sinh(2|h)cosh(2£h)] 


(118) 
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Divide the numerator and the denominator by K(1 - v) so that 


\(r, -h) = ap 


/•OO 

f 

0 (|r)Jj(a|)[sinh(24h) + 2£h cosh(2i-h)]d£ 

1 ^ 

[(2£h)2 - sinh^(2£hj] 

- [2ijh + sinh(2|h)cosh(2£h)] 

0 

£ 

H* 

J 


(119) 


When the following substitutions p = r/2h, a = a/2h, and 77 = 2h| are made, equa- 
tion (115) may be written as 




J q(p t ?) J j( a ??) [sinh(T7 ) + 77 cosh(7j)]d77 
I 2hK(l ?? - v) & 2 " sinh2 fo3 “ N + sinh(77)cosh(?7)] 


( 120 ) 


Let G = 


2hK(l - -yj and p = F/7ra“; therefore, 


= F/4a 2 ; 


a z (r - - h > = T- 5 


-/ 


J 0 (P T ?)J 1 («T7)[sinh(7j) + 77COsh(77)]d77 
477 h^a^Q Gt?^ - sinh^^] - [77 + sinh( 77 )cosh( 7 ])] 


( 121 ) 


a nh 

The dimensional quantity l — =— I has been plotted in figure 23 for various 

' r, -h 

values of a, p, and G. 


Discussion of Analytical Results 


The curves in figure 23 show the distribution of contact stress a at z = -h for 

z 


various foundation moduli and circle radii. Curves for a 


a , u, and w 
rz’ ’ 


z’ r’ 

could have been plotted for any I z| < h. Because w = cr^K at the contact surface, 

the curves also represent the deflection in the z-direction and, therefore, permit the 
visualization of the deflection at that surface. 


Because equation (121) is rather complicated to evaluate, it was necessary to use 
numerical integration techniques. The equation was evaluated on the Univac 1108 com- 
puter by the use of the Gaussian integration formula. Because the limits of integration 
are zero to infinity, it is necessary to perform the integration repeatedly with limits of 
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zero to lOn for n = 1, 2, 3, , , until an increase in n would cause less than 0.00001 
increase in the accuracy when compared to the last integration. Each calculated point 
of each curve a = 0, a = h, a = 2h for each value of K is computed in this way. 

A few statements must be made about the curves in general. Note that each 
curve crosses the r/2h axis. This means that the elastic foundation is actually put 
in tension at certain places and that the displacement in the z-direction actually be- 
comes positive. This fact is not shown in the curves presented by Sneddon (ref. 21) 
for a rigid foundation. As a further check on the validity of the solution, any of the 
curves presented in this section can be used to show that the total reaction on the bot- 
tom layer is equal to the total applied load F; this must be the case for the true solu- 
tion. 


A few statements should be made about particular curves. For either a soft 
foundation or for a thin plate (fig, 23(a)), the contact stress is almost independent of 
the area of load distribution. For slightly larger foundation moduli (figs. 23(b) and 
23(c)), the contact stress is affected by the circle radius only in the near vicinity of the 
load. For higher ratios of 2h(l - v)K/il, the contact stress is significantly affected at 
all points along the radius. By distributing the surface load over a larger area, the 
contact stress decreases in intensity but also is spread over a larger area, as shown 
by the curves. As the ratio 2h(l - v)K/p. increases, these curves approach the curves 
for a layer resting on a rigid foundation (ref , 21). When K = °°, the curves correspond 
exactly to those presented by Sneddon (ref. 21) for a layer on a rigid foundation. 

When the numerical values for the curve 2h(l - v)/\i - 1000 and those for K = 00 
are studied, there appears to be only a 0. 24-percent difference. Therefore, it can be 
concluded that if the ratio 2h(l - v)YL/ \i is greater than 1000, an error of no more than 
0. 24 percent will be encountered by the use of Sneddon’s solution for a rigid foundation. 


CONCLUSION 


Elementary problems that involve contact problems in elasticity have been pre- 
sented, primarily to clarify the term ” contact problem.” In addition, equations have 
been derived for computation of the stresses and displacements in a two-dimensional 
half plane, in a three-dimensional half space, and in a layer resting on a rigid founda- 
tion. 


A method for determination of the stresses and displacements in a thick layer 
that is surface loaded and that rests on an elastic foundation has been developed in a . 
step-by-step manner. This has provided not only exposure to the procedure used to 
derive the equations but also the equations to be used for computational purposes. 

Then, these equations were used to compute values for the special case of uniform 
pressure distribution over a circle, and the resultant curves were plotted. The follow- 
ing observations were made from the curves: (1) the elastic foundation is actually 
placed in tension at various points on the contact surface, (2) for a soft foundation, the 
contact pressure is almost independent of the area of load distribution, (3) for a larger 
foundation modulus, the contact pressure is affected by the area of load distribution 
only in the vicinity of the load, and (4) for a very large foundation modulus, the contact 
stress approaches that of a layer on a rigid foundation. 
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Two topics are suggested for future work involving variations of the problem of 
a layer resting on an elastic foundation. These topics are the indentation of a layer 
with a rigid punch and the presence of a variable foundation modulus. 


Manned Spacecraft Center 

National Aeronautics and Space Administration 
Houston, Texas, February 14, 1969 
914-50-10-06-72 
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Figure 3. - Thin plate resting on an elastic foundation. 
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Figure 4. - Pressure on bottom of thin plate resting on an elastic foundation 
and subjected to a concentrated load P. 
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Figure 15. - Half space deformed by a punch. 
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Rigid foundation 


Figure 18. - Elastic layer resting on rigid foundation and subjected to a 

surface pressure. 
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Figure 19. - Elastic layer resting on rigid foundation and subjected to a pressure 

uniformly distributed over a circle. 
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NASA-S-69-159 



(a) Layer on elastic foundation subject to an axially symmetrical pressure. 
Figure 22. - Thick layer on elastic foundation. 
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(b) Cross section through axis of symmetry showing coordinate system. 

Figure 22. - Concluded. 
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Figure 23. - Continued. 
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(k) K = °°. 

Figure 23. - Concluded. 
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